arXiv:1501.06709vl [astro-ph.GA] 27 Jan 2015 


Draft version January 28, 2015 

Preprint typeset using ETgX style emulateapj v. 5/2/11 


EXTREMELY METAL-POOR GALAXIES: THE ENVIRONMENT 

M. E. Filho 1,2 ’ 3,4 ’ 5 , J. Sanchez Almeida 2 - 3 , C. Munoz-Tunon 2,3 , S. E. Nuza 6 , F. Kitaura 6 and S. Hess 6 

Draft version January 28, 2015 

ABSTRACT 

We have analyzed bibliographical observational data and theoretical predictions, in order to probe the envi¬ 
ronment in which extremely metal-poor dwarf galaxies (XMPs) reside. We have assessed the H I component 
and its relation to the optical galaxy, the cosmic web type (voids, sheets, filaments and knots), the overdensity 
parameter and analyzed the nearest galaxy neighbours. The aim is to understand the role of interactions and 
cosmological accretion flows in the XMP observational properties, particularly the triggering and feeding of 
the star formation. We find that XMPs behave similarly to Blue Compact Dwarfs; they preferably populate 
low-density environments in the local Universe: ~ 60% occupy underdense regions, and ~ 75% reside in voids 
and sheets. This is more extreme than the distribution of irregular galaxies, and in contrast to those regions pre¬ 
ferred by elliptical galaxies (knots and filaments). We further find results consistent with previous observations; 
while the environment does determine the fraction of a certain galaxy type, it does not determine the overall 
observational properties. With the exception of five documented cases (four sources with companions and 
one recent merger), XMPs do not generally show signatures of major mergers and interactions; we find only 
one XMP with a companion galaxy within a distance of 100 kpc, and the H I gas in XMPs is typically well- 
behaved, demonstrating asymmetries mostly in the outskirts. We conclude that metal-poor accretion flows may 
be driving the XMP evolution. Such cosmological accretion could explain all the major XMP observational 
properties: isolation, lack of interaction/merger signatures, asymmetric optical morphology, large amounts of 
unsettled, metal-poor H I gas, metallicity inhomogeneities, and large specific star formation. 

Subject headings: galaxies: dwarf-galaxies: formation-galaxies: evolution-galaxies: interaction-galaxies: 
starburst 


1. INTRODUCTION 

The standard model of galaxy formation and evolution in¬ 
cludes hierarchical merging/clustering, where smaller struc¬ 
tures are used as building blocks of more massive galaxies, 
and evolution through accretion from cosmological gas fila¬ 
ments in the cosmic web (e.g.. Combes 2004; Dekel et al. 
2009; Dekel, Sari & Ceverino 2009; Nuza et al. 2014a; 
Sanchez Almeida et al. 2014b). Because merger events are 
generally responsible for the disruption of gas and stellar ag¬ 
gregates, they have been traditionally associated with the trig¬ 
gering of supermassive black hole accretion at the center of 
galaxies (Active Galactic Nuclei activity; Springel et al. 2005; 
Hopkins et al. 2008; Di Matteo et al. 2008), and with the gen¬ 
eration of star formation episodes (e.g., Sanders et al. 1988; 
Robaina et al. 2009). However, the contribution of cosmolog¬ 
ical cold-accretion flows to galaxy evolution in general, and 
star formation in particular, is recently gaining ground. In¬ 
deed, cosmological accretion onto galaxies of external metal- 
poor gas is predicted by numerical simulations to be the main 
mode of galaxy and disk assembly, at high, and possibly also 
at low, redshift (Birnboim & Dekel 2003; Brooks et al. 2009). 
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The theory predicts that these cosmological filaments are the 
main process by which star formation is fed and triggered at 
all cosmological times (Brooks et al. 2009; Verbeke et al. 
2014). Cosmological gas accretion should have been more 
common at high redshift, but in the local Universe, it is hy¬ 
pothesized to intervene in the evolution of individual low- 
mass galaxies in low-density environments, via small-scale 
gas streams, facilitating the arrival of gas from large distances 
(Keres et al. 2005; Dekel & Birnboim 2006; Dekel et al. 
2009; Dekel, Sari & Ceverino 2009; Brooks et al. 2009; Cev¬ 
erino, Dekel & Bournaud 2010; Ceverino et al. 2014; Nuza et 
al. 2014a). 

There has been some observational evidence for these ac¬ 
cretion flows at high redshift, from the detection of inverse 
metallicity gradients (lower metallicity in the central star¬ 
forming regions than in the periphery of the sources; Cresci 
et al. 2010) and from high-column density H I absorbers 
(e.g., van de Voort et al. 2012). However, in the local Uni¬ 
verse, there has only been indirect evidence for these events, 
mainly from the study of star formation in dwarf galaxies. 
For example, for galaxies of the same mass, the metallicity 
was shown to decrease when star formation increases (Man- 
nucci et al. 2010; Lara-Lopez et al. 2010; Perez-Montero 
et al. 2013). It was also found that metallicities in qui¬ 
escent blue compact dwarf (BCD) galaxies are higher than 
those of star-forming BCDs (Sanchez Almeida 2008, 2009). 
More recently, observations of tadpole galaxies in the Kiso 
survey of UV-bright galaxies (Elmegreen et al. 2013) and 
extremely metal-poor dwarf star-forming galaxies (XMPs), 
have shown that the regions where star formation occurs pos¬ 
sess a lower metallicity than the underlying galaxy (Sanchez 
Almeida 2013, 2014a). From the H I data, it has been shown 
that local BCDs and XMPs are generally surrounded by pools 
and streams of H I gas, over three times larger than the op- 
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tical extension, which frequently display, in the peripheries, 
distorted and lopsided H I distributions and velocities (Be¬ 
gum & Chengalur 2003; Young et al. 2003; Thuan, Hibbard 
& Levrier 2004; Ekta, Chengalur & Pustilnik 2006; Begum 
et al. 2006; Chengalur et al. 2006; Walter et al. 2008; Ekta, 
Chengalur & Pustilnik 2008, 2009; Ekta & Chengalur 2010; 
Ott et al. 2012; Lelli et al. 2012a, b; Lelli; Verheijan & Frater- 
nali 2014a, b; Beccari et al. 2014). Furthermore, significant 
quantities of metal-poor, large-scale H I gas were found in the 
subset of XMPs (Filho et al 2013; hereinafter Paper I). Be¬ 
cause the gas in galaxy disks is diluted in a relatively small 
timescale (of the order of one galaxy rotation; e.g., de Avillez 
& Mac Low 2002; Yang & Krumholz 2012), this cumula¬ 
tive evidence suggests that the gas had to have been accreted 
recently, most likely as vestigial streams of cosmological ac¬ 
cretion (Sanchez Almeida 2014b), like the streams modeled 
in high-mass galaxies (Dekel et al 2009; Dekel, Sari & Cev- 
erino 2009; Ceverino, Dekel & Bournaud 2010; Ceverino et 
al. 2014). The rationale is that the metallicity of the accreted 
H I gas is low. As this gas is fed to the star formation regions 
of the source, the overall metallicity in this region is lowered, 
due to dilution of the gas, while elsewhere in the galaxy, the 
H I gas has been enriched with metals from stellar winds and 
supernovae from previous star formation episodes. 

Thus, XMPs may constitute the ideal local laboratories to 
probe the role of cosmological accretion flows in galaxy evo¬ 
lution and star formation. XMPs comprise over 0.1% of 
the galaxies in the local volume (Morales-Luis et al. 2011), 
and represent the low-metallicity end of the most metal-poor 
galaxies, below an arbitrarily defined metallicity of one-tenth 
solar, or 12 + log(0/H) < 7.65. Although most XMPs are 
a subsample of star-forming BCDs (Sargent & Searle 1970; 
Thuan & Martin 1981; Papaderos et al. 1996a; Papaderos 
et al. 1996b; Telles & Terlevich 1997; Kunth & Ostlin 2000; 
Cairos et al. 2001; Bergvall & Ostlin 2002; Cairos et al. 2003; 
Noeske et al. 2003; Gil de Paz & Madore 2005; Amorfn et al. 

2007, 2009; Micheva et al. 2013), some XMPs are more dif¬ 
fuse and correspond to the category of dwarf irregular galax¬ 
ies (e.g., Skillman et al. 2013). They contain copious amounts 
of metal-poor, unsettled, H I gas (Begum & Chengalur 2003; 
Young et al. 2003; Thuan, Hibbard & Levrier 2004; Ekta, 
Chengalur & Pustilnik 2006; Begum et al. 2006; Chengalur 
et al. 2006; Walter et al. 2008; Ekta, Chengalur & Pustilnik 

2008, 2009; Ekta & Chengalur 2010; Ott et al. 2012; Lelli 
et al. 2012a, b; Paper I; Lelli, Verheijan & Fraternali 2014a, 
b; Beccari et al. 2014). Moreover, in ~ 80% of the cases, 
XMPs show evidence for optical asymmetry (Papaderos et al. 
2008; Morales-Luis et al. 2011; Paper I). The lopsided optical 
morphology has been commonly interpreted as asymmetric 
star formation, occurring in a galaxy/disk that is in a transient 
phase (Papaderos et al. 2008; Elmegreen 2009; Elmegreen & 
Elmegreen 2010), a rare occurrence at low redshift. 

The objective of this paper (Paper II) is to analyze the en¬ 
vironment of the XMPs, with the aim of understanding the 
role of cosmological accretion flows in the galaxy/disk as¬ 
sembly and triggering/feeding of the star formation in XMPs. 
Archival high resolution interferometric H I data are used 
to investigate the nearby environment, whereas cosmologi¬ 
cal simulations constrained by observations, and the observed 
distribution of galaxies, provide information on the environ¬ 
ment on large scales. The sources are the same as in Paper I, 
and haven been taken from the parent sample of 140 local 
XMPs, selected by Morales-Luis et al. (2011), from the Sloan 


Digital Sky Survey (SDSS) Data Release (DR) 7 (Abazajian 
et al. 2009) and literature. For 19 XMPs, archive H I imaging, 
at resolutions higher than 1', are available via interferomet¬ 
ric Very Large Array (VLA) and/or Giant Metrewave Radio 
Telescope (GMRT) observations; these are used to probe the 
nearby environment. The cosmic web environment is inferred 
from the constrained numerical simulations by HeB, Kitaura 
& Gottlober (2013) and Nuza et al. (2014b), and the catalog 
of SDSS galaxy filaments assembled by Tempel et al. (2014). 
We have also searched for close companions to the XMPs in 
the SDSS database. 

The paper is organized as follows. In Section 2 we present 
a summary of the optical and H I properties of XMPs with 
archival high resolution H I observations, which includes an 
analysis of the small-scale environment. Section 3 contains 
our investigation of the large-scale environment of XMPs in 
light of three diagnostic tools: cosmic web type, overdensity 
and nearest neighbours. In Section 4 we discuss the results 
and present the conclusions. The Appendix contains a quali¬ 
tative description of the archival high resolution H I morphol¬ 
ogy and velocity field of individual XMPs. 

2. SMALL-SCALE ENVIRONMENT 

Of the 140 local XMPs in the parent sample, 29 sources 
with no H I information have been observed with the single¬ 
dish Effelsberg radio telescope, resulting in 11 new detections 
(Paper I). Overall, 53 XMPs possess H I data in literature; but 
of these, only 19 have been observed at high spatial resolu¬ 
tion (< 1'), with either the GMRT or the VLA. We use these 
archival H I observations to investigate the spatial and dynam¬ 
ical distribution of the H I gas, and its relation to the optical 
component, which can provide clues as to the origin of the 
observed XMP properties. 

2.1. XMPs With Archival High Resolution H I Observations 

Lelli, Verheijan & Fraternali (2014a, b) have presented a 
detailed quantitative assessment of the H I properties (from 
new and archival data) of 18 starburst dwarf galaxies, four of 
which are also in our XMP sample (UGC 4483, UGC 6456, 
IZw 18, and SBS 1415+437). In this paper, we attempt only 
a qualitative description of the high resolution H I data in 
literature, in order to obtain global statistical properties for 
the XMP sample. 

In Table 1 we present a summary of the optical properties 
of the sample of 19 local XMPs with archival high resolution 
interferometric H I observations. The data is adapted from 
Table 3 of Paper I, with some revision and description of the 
optical morphology. Table 2 provides a summary of the rel¬ 
evant H I properties, at several scales, of the 19 local XMPs 
with archival high resolution interferometric H I observations. 
The properties are a compilation of interferometric H I data 
from literature, namely VLA and GMRT data. Low resolution 
maps display resolutions > 30", while intermediate resolution 
display 30 - 15" and high resolution display < 15". In the 
Appendix, we present a detailed description of the H I high 
resolution data for each XMP. 

If we compare the XMP H I morphology and velocity field 
with the classification given in Lelli, Verheijan & Fraternali 
(2014b), we conclude that almost all of the XMPs fall within 
their classification of Type A/B sources - galaxies with regu¬ 
larly rotating H I disks/galaxies with kinematically disturbed 
H I disks. That is, the XMPs generally show an H I disk- 
type morphology, with a velocity field typical of a rotating 
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TABLE 1 

Summary of the optical data for the 19 extremely metal-poor galaxies in the local Universe with archival high 

RESOLUTION INTERFEROMETRIC HI OBSERVATIONS. 


Source 

Name 

(1) 

RA(J2000) 

h m s 

(2) 

DEC(J2000) 

o t rt 

(3) 

m g 

mag 

(4) 

^opt 

(5) 

PA 

(6) 

O/H 

(7) 

Optical 

Morphology 

(8) 

Note 

(9) 

J0119-0935 

01 19 14 

-09 35 46 

19.5 

21.8 

115 

7.31 

cometary 


UM 133 

01 44 42 

+04 53 42 

15.4 

54.7 

105 

7.63 

cometary 

1 

SBS 0335-052W 

03 37 38 

-05 02 37 

19.0 

5.0 


7.11 

symmetric 


SBS 0335-052E 

03 37 44 

-05 02 40 

16.3 

8.3 

45 

7.31 

cometary 


UGC 4305 

08 19 05 

+70 43 12 


26.5 

115 

7.65 

multi-knot disk 

2 

HS 0822+03542 

08 25 55 

+35 32 31 

17.8 

13.9 

45 

7.35 

cometary 


DD 053 

08 34 07 

+66 10 54 

20.3 

71.2 

55 

7.62 

multi-knot disk 

3 

UGC 4483 

08 37 03 

+69 46 31 

15.1 

61.8 

80 

7.58 

cometary 

4 

IZw 18 

09 34 02 

+55 14 25 

16.4 

18.5 

45 

7.17 

two-knot cometary 

5 

Leo A 

09 59 26 

+30 44 47 

19.0 

272.0 

10 

7.30 

multi-knot disk 


Sextans B 

10 00 00 

+05 19 56 

20.5 

269.0 

25 

7.53 

multi-knot disk 

6 

Sextans A 

10 11 00 

-04 41 34 


264.0 

45 

7.54 

multi-knot disk/symmetric 

7 

UGC 6456 

11 28 00 

+78 59 39 


93.0 

115 

7.35 

multi-knot disk 

8 

SBS 1129+576 

11 32 02 

+57 22 46 

16.7 

45.6 

75 

7.36 

cometary 


UGCA 292 

12 38 40 

+32 46 01 

18.9 

57.0 

165 

7.28 

multi-knot disk 

9 

GR 8 

12 58 40 

+ 14 13 03 

17.9 

78.0 

130 

7.65 

multi-knot disk 

10 

SBS 1415+437 

14 17 01 

+43 30 05 

17.8 

55.5 

120 

7.43 

two-knot cometary 

11 

Sag Dig 

19 29 59 

-17 40 41 


116.0 

170 

7.44 

multi-knot disk/symmetric 

12 

J2104-0035 

21 04 55 

-00 35 22 

17.9 

13.1 

110 

7.05 

cometary 



Columns: 

Col. 1: Source name. 

Col. 2: Right Ascension (J2000). 

Col. 3: Declination (J2000). 

Col. 4: SDSS DR7 g-band (Petrosian) magnitude, from Morales-Luis et al. (2011). 

Col. 5: Optical diameter along the major axis, obtained from the SDSS DR7 r-band or Digitized Sky Survey (DSS)-II 
7?-band data, measured at 25 mag arcsec -2 . 

Col. 6: Approximate optical position angle, along the largest extension, measured counterclockwise from the West, 
from the SDSS DR 10 or DSS images. 

Col. 7: Metallicity, 12 + log (O/H), from Morales-Luis et al. (2011). 

Col. 8: Optical morphology, obtained from the SDSS DR10 or DSS images. Symmetric (overall spherically symmetric 
structure), disk (overall disk-like structure), two-knot (structure with two prominent star formation knots), multi-knot 
(structure with multiple star formation knots) and cometary (overall head-tail structure, with at least one prominent 
star formation knot), are used. 

Col. 9: Note on the optical structure from the SDSS DR10 or DSS images. 

Note in column 9: 

1 - The extended cometary tail is curved to the West, in a ”boomerang-type” structure. 

2 - There is a faint nodule of emission to the Northwest, a bright ’’loop-shape” structure to the East, and a faint 
filament, beginning near the loop, and curving to the Southwest. The brightest star formation knots are found on the 
Eastern edge of disk, near the loop and filament. 

3 - There is second disk-like emission to the West, at a position angle of approximately 55°, similar to the main disk. 
The brightest star formation knots are found on the Northwestern tip of the main disk. 

4 - Although the overall shape is cometary, there appears to be a fainter X-shape pattern, with the brighter and the 
fainter arm at a position angle of approximately 95 and 60°, respectively. 

5 - The two brightest star formation knots are each at the end of the cometary structure. 

6 - There is a faint tail of emission to the West, at a position angle of approximately 160°. 

7 - The brightest star formation knots are found on the Eastern tip of the disk/symmetric structure. 

8 - There is faint extended emission along a position angle of approximately 75°. The brightest star formation knots 
are found on the Eastern edge of the disk. 

9 - The disk has a ’’boomerang” shape, with the brightest star formation knots in the center and on the Eastern tip of 
the disk. 

10 - The disk shows curvature at both ends, in a ”boomerang-type” morphology. 

11 - The two brightest star formation knots are at the head of the cometary structure. The extended cometary tail is 
curved to the West. 

12 - The brightest star formation knots are found on the Eastern tip of the disk/symmetric structure. 


disk, which may present some small asymmetries and irregu¬ 
larities in the outskirts, and which may be (only slightly) of- 
fcenter relative to the optical peak and/or star-forming knots; 
statistically, we do not observe significant perturbations in the 
large-scale H I gas component of the XMP population (Ap¬ 
pendix). High spatial resolution observations of the H I gas 
suggest feeding of the star formation regions; the gas is re¬ 
solved into complex substructure, with the gas ’’pointing at” 


or ’’cradling” bright star-forming regions, or showing signs of 
’’clearing” around the star-formation knots (Appendix). Sag 
Dig has been suggested to be the result of a merger process 
(Beccari et al. 2014 and Appendix). Its H I morphology 
is ring-like, with an offcenter depression, and the velocity 
field displays no sign of rotation. These features are gener¬ 
ally not found in XMPs, which we interpret as evidence that 
the XMPs, as a class, are currently not undergoing a major 



4 


M. E. Filho et al. 


TABLE 2 

Summary of the H I data of the 19 extremely metal-poor galaxies in the local Universe with archival high 

RESOLUTION INTERFEROMETRIC HI OBSERVATIONS. 


Source 

Name 

(i) 

PA 0 pt 

(2) 

HI/opt Ratio 

(3) 

HI/opt Offset 

(4) 

PAhi 

(5) 

PAvg 

(6) 

Note 

(7) 

Ref 

(8) 

J0119-935 

115 

4 

10/10 

120-110 

120 


1 

UM 133 

105 

2 

0/10 

105 

105 


1 

SBS 0335-052W“ 


6 

10/10 

150 

NE-SW-N-S 

P 

2 

SBS 0335-052E“ 

45 

4 

10/10 

145-20 

NE-SW 

P 

2 

UGC 4305 

115 

5 

0/0 


100 


3 

HS 0S22+03542 

45 

10 

10/10 

70-80 

70 

P 

4 

DD 053 

55 

3 

10/10 

150 

45 


3,5 

UGC 4483 

80 

2 

10/0 

80-60-80 

90-120 


6,7,8,9 

IZw 18 

45 

7 

10/0 

45 

50-55 

P 

8,9,10 

Leo A 

10 

3 


10 



5 

Sextans B 

25 

3 

10/10 

25 

130 


6 

Sextans A 

45 

2 

10/0 

135 

150 


6 

UGC 6456 

115 

2 

10/0 

75 

90 


8,9,11 

SBS 1129+576 

75 

2 

0/0 

75-70 

70 

P 

12 

UGCA 292 

165 

3 

0/0 

175 

145 


6, 13 

GR 8 

130 

2 

10/10 

165 

80 


5,6,13,14 

SBS 1415+437 

120 

2 

10/0 

115 

115 


8,9 

Sag Dig 

170 

6 

0/0 



M 

5, 15 

J2104-0035 

110 

5 

10/10 

90 

90 


16 


a SBS 0335-052E and SBS 0335-052W are an interacting pair. Due to lack of resolution, the PAvg is 
merely indicative. 


Columns: 

Col. 1: Source name. 

Col. 2: Approximate position angle of the largest optical extension, measured counterclockwise from the West, 
from the SDSS DR10 or DSS images. 

Col. 3: Approximate Hl-to-optical ratio, where the HI extension is measured in the lowest resolution interfero¬ 
metric map. 

Col. 4: Approximate offset between the HI peak and the optical center/brightest star-forming knot. 

Col. 5: Approximate position angle of the largest HI extension, measured counterclockwise from the West (from 
the lowest to highest resolution map). 

Col. 6: Approximate position angle of the HI velocity gradient, measured counterclockwise from the West (from 
the lowest to highest resolution map). 

Col. 7: Note on the possible interaction/merger of the XMP. ”P” designates a galaxy that belongs to a confirmed 
interacting pair. ”M” designates a galaxy that may have undergone a recent merger. 

Col. 8: Reference, with interferometer data and/or survey. 

Instrument, interferometer and survey acronym: 

GMRT - Giant Metrewave Radio Telescope 
VLA - Very Large Array 
THINGS - The HI Nearby Galaxy Survey 
ACS - Advanced Camera for Surveys 
ANGST - ACS Nearby Galaxy Survey Treasury 

Reference in column 8: 

1 - Ekta & Chengalur (2010); GMRT 

2 - Ekta, Pustilnik & Chengalur (2009); GMRT 

3 - Walter et al. (2008); VLA THINGS 

4 - Chengalur et al. (2006); GMRT 

5 - Begum et al. (2006); GMRT 

6 - Ott et al. (2012); VLA ANGST 

7 - Lelli et al. (2012b); VLA 

8 - Lelli, Verheijan & Fratemali (2014a); VLA 

9 - Lelli; Verheijan & Fratemali (2014b); VLA 

10 - Lelli et al. (2012a); VLA 

11 - Thuan, Hibbard & Levrier (2004); VLA 

12 - Ekta, Chengalur & Pustilnik (2006); GMRT 

13 - Young et al. (2003); VLA 

14 - Begum & Chengalur (2003); VLA 

15 - Beccari et al. (2014); VLA 

16 - Ekta, Chengalur & Pustilnik (2008); GMRT 
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interaction or merger. Rather, the asymmetries and irregular¬ 
ities of the outer envelope are understood as signs of weaker 
interactions caused by minor mergers, gas accretion or gas 
outflows driven by galactic winds. 

3. LARGE-SCALE ENVIRONMENT 

In order to investigate the environment in which XMPs oc¬ 
cur, and the influence that cosmological accretion flows may 
have on their evolution and observed properties, we have com¬ 
puted the likelihood of an XMP to inhabit a certain type of 
environment using three diagnostics: cosmic web type, over¬ 
density and nearest neighbours. 

The environment can be classified according to local indi¬ 
cators, such as the local density, or according to non-local 
indicators, based, for example, on the tidal-field tensor. In the 
first case, one can estimate the local overdensity based solely 
on observations, defined as the number density of galaxies 
relative to the background number density, or by characteriz¬ 
ing the nearest neighbours (type, mass, number and distance). 
However, a classification of the environment into its cosmic 
web components requires knowledge of the large-scale envi¬ 
ronment. The sign of the eigenvalues of the tidal-field tensor 
indicate pure contracting regions, such as clusters or knots, 
or expanding regions in three, two, or one dimensions, rep¬ 
resenting voids, sheets, or filaments, respectively (Hahn et al. 
2007). While numerical simulations provide all the necessary 
information to perform a cosmic web analysis, observations 
require a prior reconstruction of the large-scale matter density 
field. The particular reconstruction technique will determine 
the accuracy of the cosmic web classification. 

3.1. Cosmic Web Type and Overdensity 

We here rely on precise constrained N-body simulations of 
the Local Universe, using the Two-Micron All-Sky Galaxy 
Redshift Survey (2MRS; Huchra et al. 2012), within volumes 
of 180 li 5 Mpc 3 (Hell Kitaura & Gottlober 2013). Cosmo¬ 
logical parameters only provide a statistical description of the 
Universe, so that an infinite number of realizations are con¬ 
sistent with the same set of parameters. The procedure pro¬ 
vides a cosmological simulation, which reproduces, in detail, 
the spatial distribution of galaxies in the local universe as de¬ 
scribed by 2MRS. The method to reconstruct the initial fluctu¬ 
ations used in these simulations is based on a self-consistent 
Bayesian machine-learning algorithm, reaching an accuracy 
of ~ 2 Mpc /r 1 (Kitaura et al. 2012; Kitaura 2013; Hel.i, Ki¬ 
taura & Gottlober 2013). The cosmic web classification used 
in this study has been performed following Forero-Romero et 
al. (2009), based on the best correlated constrained N-body 
simulation with the 2MRS galaxy field (Nuza et al. 2014b). 
The overdensity (<5) used in this study has been computed 
from the same constrained simulation. The preference of a 
type of galaxy to be located in a particular environment is 
quantified by // (r, e), an excess probability ratio, where r is 
the galaxy type and e is the cosmic web environment type 
(voids, sheets, filaments and knots; Nuza et al. 2014b). The 
environment (overdensity or cosmic web type) of a source can 
be assigned according to the cell in the reconstruction where 
the source is located. Each cell is characterized by three co¬ 
ordinates: right ascension, declination, and redshift, with the 
redshift parameterizing the distance, including the proper mo¬ 
tions caused by local tidal fields, taken consistently into ac¬ 
count. 

Because we possess coordinates for our XMP galaxies, we 
have used the constrained simulation to determine the type of 


environment (comic web type and overdensity) in which the 
XMPs reside. Spectroscopic redshifts from the SDSS DR10 
(Ahn et al. 2012) are available for ~ 70% of the 140 XMPs. 
For the remaining XMPs, photometric redshifts from the 
SDSS DR10 or NASA/IPAC Extragalactic Databasefl (NED) 
were used. 

In Figure 1, we have compared the cosmic web type en¬ 
vironment distribution of XMPs with other galaxy morpho¬ 
logical types, namely, ellipticals (E), lenticulares (SO), spirals 
(Sp) and irregular galaxies (Irrs), as assigned by the 2MRS 
team (Huchra et al. 2012). Although the Poisson errorbars 
are larger for the XMP sample, it is apparent that the XMPs 
follow similar trends as Irr galaxies, but with a more pro¬ 
nounced behaviour (Fig. 1). The XMP prevalence in voids 
and sheets, and the avoidance of knots, suggests that, statisti¬ 
cally and on a large scale (~ 1 Mpc), the XMPs are relatively 
isolated (Fig. 1). 

As a consistency check, we have also estimated the cosmic 
web type distribution for a large sample of control galaxies (~ 
700 BCDs, with Z > 0.1 Z 0 ) from the SDSS DR8 (Aihara et 
al. 2011). We find a similar behaviour among the BCDs and 
XMP sample, which is markedly different from the environ¬ 
ments found for E, SO and Sp galaxies, and more extreme than 
the behaviour of the Irr galaxies in the zero redshift 2MRS 
sample (to be included in a future paper). 

Figure 2 shows the projected spatial distribution of the 140 
local XMPs, color-coded for cosmic web type. The sizes 
of the circles trace the overdensity (<5), such that, the larger 
the circle, the larger the overdensity. The values range from 
the most underdense (5 = -0.9, for PHL 293B) to the most 
overdense (<5 = 376, for VCC 0428, a BCD in the Virgo clus¬ 
ter; Meyer et al. 2014) source environment. Figure 3 shows 
the overdensity distribution (log(l+<5)), color-coded by cos¬ 
mic web type, for the XMP sample. Arrows, color-coded by 
cosmic web type, show the mean overdensity (log( 1 + < 6 >)) 
for each distribution. According to the overdensity distribu¬ 
tion, ~ 60% of the XMPs are found in underdense regions 
(5 < 0) of the reconstruction (Fig. 2 and 3). On the other 
hand, most of the XMP galaxies (~ 75%) reside in the less 
dense cosmic web type environments (~ 40% are in sheets 
and ~ 34% in voids), while ~ 25% can be found in filaments, 
and < 1% in knots (Fig. 1, 2 and 3). 

In the analysis below, because we are constrained by small 
number statistics and systematic errors that may affect the de¬ 
termination of the physical parameters, the trends should be 
taken with caution. 

We have investigated whether the H I mass (Mhi), dynam¬ 
ical mass (M Dyn ) or stellar mass (M.) of the XMPs depends 
on the environment in which they reside (Fig. 4). The H I 
and stellar masses were taken from Table 5 of Paper I. In this 
Paper II, we have estimated the dynamical masses with re¬ 
vised optical radius measurements. There is a large overlap 
in dynamical, H I and stellar mass for a range of overdensi¬ 
ties (Fig. 4). However, XMPs in filaments and sheets show 
the largest range in H I mass, particularly extending to lower 
H 1 masses, while XMPs in voids show tendentiously larger 
H I masses, above log (Mhi/M q ) =* 8 (Fig. 4a). Five sources 
in voids and sheets show the largest dynamical masses (log 
(Moyn/Mo) > 9; Fig. 4b). In terms of stellar mass, most XMPs 
in filaments found in overdense regions (6 > 0) show stellar 
masses below log (M*/M Q ) - 7, with a tail of XMPs in less 
dense filamentary environments towards high stellar masses 

7 http://ned.ipac.caltech.edu/ 
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Fig. 1.— The excess probability distribution of XMPs, and E, SO, Sp and Irr galaxies of the zero redshift 2MRS sample (Huchra et al. 2012; Nuza et al. 2014b). 
The four panels correspond to sources found in voids, sheets, filaments and knots. 
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Fig. 2.— The projected spatial distribution of the 140 local XMPs, color- 
coded by cosmic web type (blue for voids, red for sheets, green for filaments 
and cyan for knots). The sizes of the circles trace the overdensity, such that, 
the larger the circle, the larger the overdensity. The values range from the 
most underdense (6 = -0.9, for PHL 293B) to the most overdense (6 = 376, 
for VCC 0428, a BCD in the Virgo cluster) source environment. 


Fig. 3.— The overdensity distribution for the 140 local XMP sample, color- 
coded by cosmic web type (blue for voids, red for sheets, green for filaments 
and cyan for knots). Arrows, color-coded by cosmic web type, show the mean 
overdensity (log(l + < 6 >)) for each distribution. 


(log (M*/M 0 ) > 8; Fig. 4c). The largest H I-to-stellar mass 
ratio belongs to an XMP residing in a void (Fig. 4d). Overall, 
we see no clear relation between the environment in which 
XMPs reside and their dynamical, stellar or H I mass. 

We have further investigated if there is an environment 
trend with metallicity (12 + log(0/H)), star formation rate 
(SFR), specific star formation rate (sSFR = SFR/M*), H I gas 
consumption rate (l/tM H i = SFR/Mhi), and optical morphol¬ 
ogy (symmetric, cometary, two-knot and multi-knot; Fig. 5). 
These values were taken from Table 3 and 4 of Paper I. We 
find no clear correlation between environment and metallicity 


(Fig. 5a), SFR (Fig. 5b) and gas consumption rate (Fig. 5d). 
However, XMPs in voids show a slight tendency for larger 
values of sSFR (log sSFR > -8 yr _1 ; Fig. 5c). Cometary and 
multi-knot XMPs are roughly equally (fractionally) prevalent 
in voids, sheets and filaments (Fig. 5e and 5f), while two-knot 
sources are roughly equally (fractionally) prevalent in voids 
and sheets (Fig. 5e and 5f). Symmetric sources are found to 
be more common in voids (Fig. 5e and 5f). The two sources 
with morphological information, and contained in knots, are 
both cometary (Fig. 5d and 5e). 

In summary, we find that XMPs tend to reside in low- 
density environments, and in cosmic web types characterized 
as voids and sheets. We further find that XMPs, BCDs and Irr 
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galaxies tend to inhabit the same type of environment, which 
is in stark contrast with the behaviour of E, SO, Sp galaxies in 
the zero redshift 2MRS sample. 

3.2. Proximity to SDSS Morphological Filaments 

We have investigated the proximity (or inclusion) of XMPs 
to morphological filaments, i.e., ’’chains” of galaxies, as em¬ 
pirically defined using galaxies in the spectroscopic sample of 
the SDSS DR8 dataset. Tempel et al. (2014) use a marked- 
point process with interactions, called the Bisous model (Sto- 
ica, Gregori & Mateu 2005), to detect galaxy filaments in the 
spectroscopic sample of the SDSS DR8 dataset (Tempel, Tago 
& Liivamagi 2012). A redshift interval of 0.009 < z < 0.155 
was used so that the lower redshift limit excludes the Local 
Supercluster. The authors look for morphological filaments 
with a radius of 0.5 h~ l Mpc. They find that the shortest fil¬ 
aments in their filament catalog are 3 — 5 /? 1 Mpc in length, 
while typical filament lengths are of 60 Mpc. 

We have mined the catalog of galaxies (hereinafter A3), 
used to generate the morphological filaments, to search for 
XMP identifications, or, alternatively, to find the nearest 
galaxy to each XMP. Distances between the XMPs and A3 
galaxies were estimated neglecting galaxy proper motions, 
since the local flow for pairs of nearby galaxies should be of 
the same order. A3 also includes the distance of each galaxy 
to the nearest morphological filament (Tempel et al. 2014). 


Hence, the distance of the nearest galaxy and morphological 
filament can be used as a proxy for the distance of the XMPs 
to morphological filaments. In order to identify our sample 
XMPs in the A3, we searched the A3 for galaxies within D ~ 
10 kpc of the XMP positions; none were found. We are, there¬ 
fore, confident that the nearest galaxy identified for each XMP 
is truly a separate galaxy. We stress that A3 contains ~ 500 
000 SDSS DR8 sources in a contiguous field, which could 
miss the XMPs and some of their possible neighbours (see 
below). 

Figure 6 contains the XMP distance to the nearest galaxy 
(Dg a iaxy) in A3, as a function of the distance between the 
nearest galaxy in A3 and the nearest morphological filament 
(^filament), color-coded by cosmic web type of the XMP. The 
figure also contains solid black lines, indicating the radii 
of the morphological filaments, so that XMPs in galaxy fil¬ 
aments should display D ga i axy < 0.5 Mpc and Dfji ament < 
0.5 Mpc. Only three of the XMPs are found within the 
f^gaiaxy 5; 0.5 Mpc and /-^lament £ 0.5 Mpc region, and only 
marginally (Fig. 6). Therefore, we find that XMPs are gener¬ 
ally not contained within the morphological filaments defined 
by Tempel et al. (2014). 

3.3. Nearest Galaxies and Potential Perturhers 

According to the galaxies in the A3 catalog, only two XMPs 
have a galaxy within D ga i axy <100 kpc, with one of these 
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Fig. 5.— The overdensity as a function of (a) metallicity, (b) SFR, (c) sSFR, (d) H I gas consumption rate (SFR/Mhi), and (e and f) optical morphology 
(symmetric, cometary, two-knot and multi-knot), color-coded by cosmic web type (blue for voids, red for sheets, green for filaments and cyan for knots). These 
parameters were taken from Table 3 and 4 of Paper I. 


XMPs residing in a filament and another in a void (Fig. 6). 
Indeed, there are four documented XMPs that are known to be 
interacting with nearby companions (IZw 18, SBS 0335-052, 
HS 0822+03542 and SBS 1129+576; Sect. 2 and Appendix). 

Overall, the nearest galaxies in A3 show a large range in 
distances to the XMPs (Z) ga iaxy ~ 0.1 - 100 Mpc), which 
does not seem to depend strongly on the XMP environment 
(Fig. 6). However, there is a tendency for XMPs in filaments 
to show a somewhat bimodal behaviour: there is a group of 
(nearest) galaxies at smaller (D ga i axy < 1 Mpc) and larger 


(^galaxy ~ 10 Mpc) distances (Fig. 6). 

In order to complement the results based on A3, we have 
also looked for nearest galaxy neighbours, as potential per- 
turbers, using the full SDSS DR9 (Ahn et al. 2012) catalog. 
The search was performed within a 12' radius around each 
XMP, corresponding to a distance of D ~ 200 kpc at z = 0.01. 
We further constrained the galaxies to have spectroscopic red- 
shifts of z <0.1, resulting in ~ 1 000 galaxies, potential per- 
turbers to 104 XMPs. We have also estimated the stellar mass 
of the potential perturbers from their SDSS DR9 g-r colours. 
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Fig. 6.— The distance between the XMP and the nearest galaxy in A3, as 
a function of the distance between the nearest galaxy in A3 and the nearest 
morphological filament, color-coded by cosmic web type of the XMP. The 
solid black lines indicate the radii of the morphological filaments, so that 
XMPs in galaxy filaments should display Z) ga j aX y 0.5 Mpc and filament ^ 
0.5 Mpc. The black dashed line indicates a representative distance for the 
nearest galaxy, in order to exert some gravitational influence over the XMP 
(see Sect. 3.3 and 3.4). 


Fig. 7.— The distance to the nearest potential perturber in the SDSS DR9 
(within D < 300 Mpc of the XMPs), as a function of stellar mass, color-coded 
by the cosmic web type of the XMP (blue for voids, red for sheets, green 
for filaments and cyan for knots). The stellar mass of the nearest potential 
perturber has been estimated from the SDSS DR9 g — r colours, using the 
mass-to-light ratios from Bell & de Jong (2001). 


using the mass-to-light ratios from Bell & de Jong (2001). 

In Figure 7 we plot the stellar mass (M*) and distance 
(Uperturber) of the nearest potential perturber within D < 
300 Mpc of the XMPs, color-coded by cosmic web type of 
the XMP Only one XMP shows a potential perturber within 
Uperturber < 100 kpc, with a stellar mass of M* ~ 10 7 M 0 
(Fig. 7). We can not, however, exclude the presence of poten¬ 
tial perturbers with no measured spectroscopic redshift in the 
SDSS, nor less luminous (and less massive) than the XMPs, 
due to the magnitude cuttoff of the SDSS (magnitude in r > 
17.7 mag). Comparing Fig. 6 and Fig. 7, we conclude that the 
potential biases in the A3 sample do not modify the conclu¬ 
sion that XMPs generally do not have close companion galax¬ 
ies. Above a D per turber ~ 1 Mpc distance to the XMP, the stel¬ 
lar masses of the nearest potential perturbers show a change; 
the stellar masses not only span a larger range, but potential 
perturbers with stellar masses of up to M« ~ 10 11 M 0 can be 
found (Fig. 7). The nearest, least massive, and most massive 
potential perturbers belong to XMPs in all three major cos¬ 
mic web types (voids, sheets and filaments). We further note 
a break in perturber mass distribution at D pemil t er ~ 10 Mpc; 
above this distance, stellar masses are skewed towards the 
larger range (M„ > 10 9 M 0 ), consistent with the detection, at 
larger distances, of more luminous, and hence, more massive 
galaxies (Fig. 7). 

Overall, with the exception of the already documented 
cases of XMPs with companions (Sect. 2 and Appendix), sta¬ 
tistically we do not find evidence for the presence of poten¬ 
tial perturbers (with measured spectroscopic redshifts in the 
SDSS) within D < 100 kpc of the XMPs. We caution, how¬ 
ever, that due to a lack of spectroscopic redshift and the spec¬ 
troscopic flux limit of the SDSS, we may be missing potential 
perturbers of very low mass and luminosity. 

3.4. Perturbations and Interactions 

Because the nearest galaxies to XMPs are generally found 
to be at distances D > 100 kpc (Fig. 6 and 7), we have at¬ 
tempted to evaluate the gravitational influence these could 
have on the XMPs, particularly if they could potentially trig¬ 


ger and/or feed the star formation. As a first approximation, 
we have analyzed the conditions for perturbation in a system 
of a satellite galaxy (XMP) on a circular orbit around a per¬ 
turber. Let us consider that the satellite galaxy has mass m and 
radius r, orbiting a perturber galaxy of mass M, at a distance 
of Z), and where m « M. 

Minor mergers/interactions (when m « M) occur through 
dynamical friction, when the perturber and satellite share a 
common dark matter halo. This dynamical friction causes 
orbital decay of the satellite galaxy within the so-called dy¬ 
namical friction timescale (e.g., Patton et al. 2000; Binney & 
Tremaine 2008): 

r fric - 2.64 • 10 5 Gyr x (1) 

min A 

where D is the distance between the satellite galaxy and the 
perturber (in kpc), v is the circular speed of the satellite (in 
km s 1 j, m is the mass of the satellite (in M 0 ), and In A is the 
Coulomb logarithm. This expression does not take into con¬ 
sideration any mass loss (from tidal stripping) that the satellite 
may suffer during the orbital decay. 

If we assume a typical distance to the perturber of D ^ 
100 kpc (Fig. 6 and 7), a satellite mass of m - 10 8 M 0 
(Fig. 4c), a Coulomb logarithm of 2 (e.g., Binney & Tremaine 
2008; Filho et al. 2014), and a circular velocity of the satel¬ 
lite of v' - sjG MID, where G is the gravitational constant, 
this then yields a dynamical friction timescale of 7j nc ^ 
50 Gyr, for M ^ m, or 150 Gyr, for M - 10 m. This esti¬ 
mate provides a timescale for the satellite to spiral in towards 
the perturber. The resulting timescales are too long for the 
merger/interaction process to be significant in XMPs. 

As the satellite approaches the perturber on this orbital de¬ 
cay, the tidal forces at work on the satellite become more im¬ 
portant. If we now consider a test particle at a distance A r 
from the surface of the satellite (where A r « r), then the 
perturber will exert a perturbative tidal force on the test parti¬ 
cle (relative to the force exerted at the center of the satellite), 
that is given by (e.g., Binney & Tremaine 2008): 
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Fig. 8.— The displacement produced by tidal forces, as a function of the 
ratio between the XMP radius and the distance to the perturber, color-coded 
by the relative XMP/perturber mass: M/m = 0.1 (cyan), 1 (blue), 10 (red), 
100 (green). The black line corresponds to A r/r = 1, the lower limit criteria 
for tidal stripping. 
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The satellite will also exert a restoring force on the test parti¬ 
cle, if tidal forces displace it from its equilibrium configura¬ 
tion at r, namely, (e.g., Binney & Tremaine 2008): 


G m A r 



(3) 


where we have neglected the pressure support in the satellite 
(Adams et al. 2011). In equilibrium Em — E m ; therefore, the 
displacement produced by tidal forces is given by: 


A r 2 M t r \ 3 
r m \d) 


(4) 


The condition for tidal disruption of the test particle by the 
perturber is producing a large perturbation, i.e., A r/r > 1. 

In Figure 8 we show the displacement produced by tidal 
forces (Eq. 4), as a function of the ratio between the XMP 
radius and the distance to the perturber, color-coded by the 
relative XMP/perturber mass. As can be seen from Figure 8, 
in order for tidal stripping to occur, the perturber is required to 
have a high mass and/or be located at a small distance to the 
satellite. For a perturber and satellite of similar mass (M - 
m ), the distance must be of the order of the size of the galaxy 
(Z) =! r). Because XMPs have radii of a few kpc, D must 
also be a few kpc in order to exert appreciable gravitational 
influence on the XMP. 

Therefore, although a significant number of XMPs show 
potential perturbers within D ~ 300 Mpc, it is unlikely that 
these will cause major gravitational disturbance (or tidal strip¬ 
ping) of the XMP, which could instigate gas transfer and sus¬ 
tain star formation. 


4. DISCUSSION AND CONCFUSIONS 

XMPs are found to tendentiously reside in low-density en¬ 
vironments (~ 60%) and in cosmic web types characterized 
as voids and sheets (~ 75%; Sect. 3.1). In this respect, their 
behaviour is extreme, in contrast to the location in the cosmic 


web of E, SO and Sp galaxies, but inhabiting regions similar 
to those of Irrs and BCDs (Sect. 3.1). 

We find results similar to those found for nearby (rela¬ 
tively) luminous galaxies (Blanton & Moustakas 2009) and 
Hcr-emitting galaxies (Darvish et al. 2014), whereby the en¬ 
vironment does not determine the overall observed properties 
(mass, luminosity, SFR, etc.) of a certain galaxy type. How¬ 
ever, the fraction of a certain galaxy type (in this case, XMPs), 
does appear to be determined by the environment in which it 
resides (Sect. 3.1). 

The evidence for occupation of low-density environments 
is further strengthened when we investigate the proximity to 
morphological filaments (Sect. 3.2), and galaxies in the vicin¬ 
ity of the XMPs (Sect. 3.3). We find no significant evidence 
for the presence of galaxy filaments within D ~ 0.5 Mpc of an 
XMP, and only one potential perturber (of mass and luminos¬ 
ity at least equal to that of an XMP, and with a measured spec¬ 
troscopic redshift in the SDSS) is found within D ~ 100 kpc 
of an XMP. However, even at these distances, we have demon¬ 
strated that a potential perturber is very unlikely to cause grav¬ 
itational disturbance capable of drawing fresh gas or driving 
the star formation (Sect. 3.4). 

We have also analyzed the small-scale environment, prob¬ 
ing for signatures of galaxy merger/interactions in the pub¬ 
lished H I data of XMPs (Sect. 2 and Appendix). IZw 18 
(Felli et al. 2012a; Felli, Verheijan & Fraternali 2014a, 
b), SBS 0335-052W and SBS 0335-052E (Ekta, Pustilnik & 
Chengalur 2009), HS 0822+03542 (Chengalur et al. 2006) 
and SBS 1129+576 (Ekta, Chengalur & Pustilnik 2006) have 
been confirmed to be interacting with companion galaxies, all 
within D <25 kpc (Sect. 2 and Appendix). Moreover, it has 
been suggested that the H I properties observed in Sag Dig 
could be the result of a recent merger between dwarf galaxies 
(Beccari et al. 2014). However, if Sag Dig is the result of a 
merger (Beccari et al. 2014; Sect. 2 and Appendix), then the 
merger signatures present in Sag Dig (ring-like H I distribu¬ 
tion versus disk-like optical morphology, clumpy H I distri¬ 
bution with no optical correspondence, large offset between 
the optical center and the depression in the H I distribution, 
no clear rotational motion) are generally absent in the XMP 
population. Indeed, we find no statistical evidence for (major) 
interactions with (gas-rich) galaxy companions in the recent 
past (Sect. 2 and Appendix; see also Brosch, Almoznino & 
Heller 2004, Holwerda et al. 2013 and Kreckel et al. 2015), a 
mechanism which is traditionally invoked as a driver for star 
formation and overall galaxy evolution (e.g.. Keel et al. 1985; 
Robaina et al. 2009). Although such an interaction could pro¬ 
vide large quantities of fresh H I gas, and create conditions for 
triggering/feeding the observed star formation, it can not ex¬ 
plain why the gas we observe is metal-poor (Paper I), nor can 
it explain the metallicity inhomogeneities observed in some 
XMPs (Sanchez Almeida et al. 2013, 2014a), and it does not 
explain why we do not systematically observe signatures of 
major disturbances/interactions in the H I gas, and particu¬ 
larly in the H I-to-optical morphological/kinematical relation 
(Sect. 2 and Appendix). In fact, the central regions of the 
HI structures are relatively well-behaved; in many cases they 
are organized in a disk with clear rotational motion, in which 
the kinematical H I axis is virtually aligned (within 15°) with 
the H I and optical axis, where the peak of the H I emission 
is positionally coincident (within 10") with the center of the 
optical disk or the brightest star formation knot (Appendix), 
and where the velocity offsets between the stellar and H I 
components are absent or small (< 40 km s' 1 ; of the order 
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of the velocity dispersion; Paper I). Hints of kinematical and 
morphological disturbance of the H I gas appear primarily in 
the galaxy outskirts; tails, extensions, clumps of H I gas, and 
kinks, discontinuities, twisting of the velocity field (Sect. 2 
and Appendix). 

Indeed using three different diagnostic tools (overdensity, 
cosmic web type and nearest galaxies/potential perturbers), 
we find evidence that XMPs are relatively isolated in the 
nearby Universe, on scales of hundreds of kpc (Sect. 2 and 
3). However, despite this isolation, we must understand how 
XMPs acquire such large amounts of metal-poor H I gas (Pa¬ 
per I and references therein), and how they sustain their star 
formation. 

An alternative scenario is that XMPs are truly primordial 
galaxies, in the sense that they have remained quiescent for 
some time, thus preserving the H I gas acquired early on in 
their formation. However, this scenario can not explain the 
current starburst phase of the XMPs, nor can it explain the 
reason, or provide a mechanism, for the XMPs to have aban¬ 
doned their quiescent phase. Moreover, because most XMPs 
have huge SFRs for their stellar masses (Fig. 5c), they would 
have had to form all their stars in less than a Gyr, at the present 
rate. 

Following Sanchez-Almeida et al. (2013, 2014b), we find 
that the XMP properties, particularly their environment char¬ 
acteristics, are consistent with XMP evolution being mainly 
driven by metal-poor gas accretion, similar to the cold accre¬ 
tion flows simulated at high redshift (Keres et al. 2005, 2009; 
Dekel et al. 2009; Dekel, Sari & Ceverino 2009; Brooks et 
al. 2009; Ceverino, Dekel & Bournaud 2010; Ceverino et al. 
2014; Nuza et al. 2014a), and supported by observational ev¬ 
idence (Cresci et al. 2010; van de Voort et al. 2012). Such ac¬ 
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high density cold gas, which are capable of penetrating low- 
mass haloes (Keres et al. 2005, 2009; Dekel et al. 2009; 
Dekel, Sari & Ceverino 2009; Brooks et al. 2009; Verbeke 
et al. 2014). Moreover, this cold-mode of accretion is able 
to explain the major observational properties of the XMPs; 
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APPENDIX 

DESCRIPTION OF XMPS WITH ARCHIVAL HIGH RESOLUTION H I OBSERVATIONS 

This appendix provides a detailed description of the H I properties of the 19 XMPs with high resolution interferometric data in 
literature. The average properties of these galaxies, as a class, are discussed in Section 2. Here we describe the peculiarities of 
the individual sources used to distill a global picture. We cite the original references of the archival data used in our analysis, as 
well as the main observational properties of the targets inferred from the visual inspection of the published maps. 

Our description of the high resolution H I properties of the XMPs are based on the definition of several source parameters: the 
position angle of the largest H I extension, measured counterclockwise from the West (PAhi), the position angle of the largest 
optical extension, measured counterclockwise from the West, from the SDSS or DSS images (PA opt ), the position angle of the H I 
velocity gradient, measured counterclockwise from the West (PAvg), the position angle of the H I velocity field (PAy n ), the offset 
between the H I peak and/or the optical center/brightest star-forming knot (AD), and the offset between the H I velocity gradient 
or largest extension position angle and the optical largest extension position angle (A0). References and survey/interferometer 
information (GMRT - Giant Metrewave Telescope, VLA - Very Large Array, THINGS - The H I Nearby Galaxy Survey, ACS 
- Advanced Camera for Surveys, ANGST - ACS Nearby Galaxy Survey Treasury) for the H I data are provided for each of the 
sources. 

J0119-935 From Ekta & Chengalur (2010; their Fig. 3 and 4), with the GMRT. The lower resolution (46" x 32") H I map 
shows a single-peaked ”boomerang-type” Northeast-Southwest disk structure (PAhi ~ 120°), with an extension to the South and 
to the North, the latter which connects to a clump of emission to the Northwest. The orientation of the H I disk is roughly 
coincident with the position angle of the optical cometary emission (PA opt » 115°). However, the center of the H I emission does 
not coincide with either the center of the optical disk, nor the bright star formation knot that leads the cometary structure in the 
Northeast (AD < 10"). In the central regions, the H I velocity gradient (PAvg ~ 120°) follows the H I and optical morphology. 
However, in the outskirts, the H I velocity field displays distortions, with the mimicking of the ’’boomerang” shape observed in 
the H I morphology. The intermediate resolution (30" x 20") image shows a ”Y-shaped” H I morphology, cradling the bright 
optical star formation knot to the Northeast, with the Southern extension curving towards the East. The position angle at this 
resolution (PAhi ~ 110°) is roughly coincident with the larger-scale H I and optical morphology. 

UM 133 From Ekta & Chengalur (2010; their Fig. 1 and 2), with the GMRT. The lower resolution (42" x 36") map of the H I 
gas displays a single-peaked Northeast-Southwest disk (PAhi ~ 105°), with a small extension to the South and a larger extension 
to the East. The H I and optical coincide in terms of position angle and geometric (disk) center/H 1 peak. There is a less than 
10" distance between the H I peak and the bright star formation knot to the Southwest, leading the cometary structure. In the 
central regions, the H I gas velocity gradient (PAvg ~ 105°) is smooth and parallel to both the optical (PA opt ~ 105°) and H I 
morphology, which is consistent with a rotating disc. However, distortions are seen in the outskirts, particularly towards the 
South and East, in the direction of the H I extensions. At an intermediate resolution (22" x 17"), the H I morphology is ’’arrow- 
shaped”, pointing towards the Northeast, and accompanying the overall optical structure. Similarly to the lower resolution 
images, extensions to the South and East are observed. The morphology and central velocity gradient on these scales, match 
the larger-scale H I and optical morphology. Although in the central regions the velocity field is well-behaved, in the periphery 
there are clear asymmetries; the velocity gradient changes to the North-Northeast, and large irregularities are detected to the 
South-Southeast. 

SBS 0335-052E and SBS 0335-052W From Ekta, Chengalur & Pustilnik (2009; their Fig. 1, 4 and 5), with the GMRT. The 
low resolution (~ 40") map shows this interacting pair of XMP galaxies (AD ~ 20 kpc), enveloped in a common H I structure 
(PAhi ~ 175°), with each H I peak corresponding to a galaxy. At intermediate (~ 20") resolution, the individual single-peaked 
H I structures are distinguishable. The elongated disks are oriented along approximately the same direction (PAhi ~ 145° for the 
Eastern, and PAhi ~ 150° for the Western, component). SBS 0335-052W displays a small nodule of emission to the Southwest. 
The H I peaks are slightly offcenter (AD < 10") from the peaks of optical emission. In the case of SBS 0335-052E, the extended 
H I emission is approximately perpendicular to the Northwestern extension seen in the optical image (PA opt « 45°). The velocity 
gradient in SBS 0035-052W changes from the center (approximately North-South) to the outskirts (Northeast-Southwest). For 
SBS 0335-052E, the velocity gradient is along the tidal tail (Northeast-Southwest), but deviates from a standard disk. High 
resolution (~ 9") observations show that SBS 0335-052W maintains the overall H I position angle seen at lower resolutions, with 
an extension to the Northeast. The direction of the velocity gradient is approximately Northwest-Southeast. SBS 0335-052E 
reveals extended emission (PAhi ~ 20°), roughly aligned with the optical cometary tail. 

UGC 4305 (also known as Holmberg II, Arp 268, DDO 050 or VII Zw 223) From Walter et al. (2008; their Fig. 13), from 
the VLA THINGS. High resolution (6.95" x 6.05") H 1 maps show an almost symmetric spiral-like structure, the internal H I 
morphology is complex, displaying holes, clumps and arms. The optical loop appears associated with the brighter H I spiral arm 
in the South, which encloses a prominent H I hole. The Southern (fainter) arm extends as a tail to the Southwest. The H I velocity 
gradient runs approximately North-South (PAvg ~ 100°), within Ad ~ 15° of the optical position angle (PA opt « 115°). 
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HS 0822+03542 From Chengalur et al. (2006; their Fig. 1, 4 and 5), with the GMRT. The low resolution (42" x 37") map 
shows the XMP target along with its pair to the Northeast, the low surface brightness galaxy SAO 0822+545, at a distance of 
D ~ 11 kpc and a position angle of PA « 125°. At this resolution, the H I gas is single-peaked, with an extension to the East- 
Northeast. The velocity field (27" x 23") displays a distortion, with the velocity gradient Northwest-Southeast (PAvg ~ 70°) near 
the center, aligned with the H I morphology (PA HI ~ 70°), and approximately East-West (PA V g ~ 170°), roughly perpendicular 
to the H I, towards the Southeast. Intermediate and high resolution (17" x 15"and 13" x 11") images show a single-peaked, 
Northwest-Southeast disk (PAhi ~ 70°), with an extension to the Northwest-West, and a clump of emission to the Northeast, the 
latter tracing the extension seen at lower resolution. The extension to the Northwest-West appears to mimic the position angle of 
the optical cometary tail (PA opt ~ 45°). There is an offset between the H I peak and the center of the optical disk and bright star 
formation knot to the Southeast (AD < 10 arcsec). At the highest resolution (7" x 5.5"), the H I gas morphology resembles two 
almost vertical (PAhi ~ 80°) "fingers”, pointing to the bright star formation knot, and peaking on either side. Overall, there is an 
offset of Ad ~ 30° between the optical and H I alignment. 

UGC 4459 (also known as DDO 053 or VII Zw 238) From Walter et al. (2008; their Fig. 17), from VLA THINGS, and 
Begum et al. (2006; their Fig. 4), with the GMRT. The H I map at intermediate resolution (29" x 27") shows a double-peaked 
quasi-symmetric spherical structure, with a small extension to the Northeast. The two H I peaks are roughly coincident with 
the opposite ends of the optical main disk, with an H I depression in between. There is slight H I enhancement associated 
with the secondary optical disk. The overall position angle of the H I emission (PAhi ~ 150°) is approximately perpendicular 
to the optical emission (PA opt ~ 55°). The velocity gradient direction (PAvg ~ 45°) is relatively well-aligned with the optical 
morphology. The velocity gradient is not uniform, with the field displaying an asymmetry, weighted towards the Southeast. 
Similarly to the intermediate resolution images, the high resolution (6.34" x 5.67") H I maps show two high-density H I clumps, 
with a central depressed region, embedded in an H I halo, with faint Eastern-Northeastern emission. The H I emission associated 
with the fainter optical disk is more apparent. On these scales, the velocity gradient direction is roughly aligned with the optical 
morphology and the larger-scale velocity gradient, and mimics what occurs at low resolution: there are deviations towards the 
Southeast. 

UGC 4483 (also known as Mrk 0090 or SBS 0826+528) From Ott et al. (2012; their Fig. 7), from the VFA ANGST, and 
Lelli et al. (2012b; their Fig. 1 and 2) and Lelli, Verheijan & Fraternali (2014a, b; their Fig. 1, C.l 1, and Fig. 1, 2), with the VLA. 
The intermediate resolution (20" x 20") H I map shows almost vertical (PAhi ~ 80°) ’’arrow-shaped” emission, where the head 
of the arrow spatially coincides with the bright star formation knot to the North, heading the cometary optical morphology. There 
is also extended H I emission to the South and Northwest. The central H I morphology matches well the observed ”X-shaped” 
and overall optical morphology (PA opt ~ 80°). The kinematical axis is vertical (FA km « 90°), roughly matching the position 
angle of the brighter arm (PA opt « 95°) associated with the ”X-shaped” optical morphology. The velocity field can be modeled 
by a rotating disk, possibly with a radial component. High resolution (10" x 10") H I observations mimic the larger-scale arrow- 
shape emission, but two H I peaks are resolved, the brightest of which is elongated to the Southeast and spatially coincides with 
the bright optical star formation knot. There is extended emission also to the South and East. On these scales, the H I major 
axis (PAhi ~ 60°) appears to follow in the direction of the fainter ’’arm” associated with the ”X-shaped” optical morphology 
(PA opt ~ 60°). The H I kinematical axis is vertical (PA kln « 90°), and matched to the larger-scale velocity gradient. The highest 
resolution (~ 6") H I image is similar to its larger-scale counterparts: the overall morphology is arrow-shaped (PAhi ~ 80°), with 
two resolved H I peaks, the brightest of which is positionally coincident with the bright star-forming knot. On these scales, the 
velocity gradient (PAvg ~ 120°) is misaligned relative to the H I and optical (PA opt « 80°) major axes. The velocity field displays 
asymmetries and distortions along the axis. 

IZw 18 (also known as UGCA 166, Mrk 0116 or SBS 0930+554) From Lelli et al. (2012a; their Fig. 2, 3 and 4) and Lelli, 
Verheijan & Fraternali (2014a, b; their Fig. C.16, and Fig. 1, 2), with the VLA. The intermediate resolution (20" x 20") H I 
image shows single-peaked ’’triangle-shaped” emission (PAhi ~ 45°), with a Southeast extension that curves due South, along 
D ~ 14 kpc. The main body of the H I emission is aligned with the optical morphology (PA opt ~ 45°). The peak of the H I 
coincides with the position of the brightest (Northwest) star formation knot in the optical image. There is small extended emission 
to the Northeast, North and Northwest. The latter emission is likely associated with the companion galaxy, IZw 18C, with which 
it is interacting, at a distance of D ~ 18 kpc. The kinematical axis iFA kln « 50°) is roughly aligned with both the (main) H I and 
optical emission. The large Southern extension is at an almost constant velocity. At high resolution (5" x 5"), the H I emission is 
resolved into two peaks, positionally coincident with the two brightest star formation knots in the optical image. There is also a 
halo of extended emission to the Northwest, englobing the interacting companion, IZw 18C. The velocity field (PA km « 55°) can 
be modeled by a rotating disk plus radial motion, with signs of distortion in the direction of IZw 18C. At even high resolution 
(2" x 2"), there is peaked H I emission extended to the South/Southeast, in an inverted ”V-shaped” morphology, and extended 
emission to the Northwest, in the direction of IZw 18C. 

Leo A (also known as Leo III, UGC 5364 or DDO 069) From Begum et al. (2006, their Fig. 9), with the VLA. The low 
resolution (~ 78") H I map shows (unequal) double-peaked H I emission, almost horizontal (PAhi ~ 10°) disk morphology, 
which is coincident with the optical position angle (PA opt « 10°). There are faint extensions to the North and South. 

Sextans B (also known as UGC 5373 or DDO 070) From Ott et al. (2012, their Fig. 10), from the VLA ANGST. The high 
resolution (~ 6") H I map shows a complex, quasi-spherical, spiral-like morphology, with H I clumps and holes. The regions 
of more enhanced H I emission appear to lie mainly above the optical disk and on the Southeast side, in an inverted ”C-shape”, 
so that the internal H I morphology is roughly aligned with the optical disk (PAhi ~ PA opt « 25°). The direction of the velocity 
gradient (PAvg ~ 130°) is virtually perpendicular to the (internal) H I and optical morphology. 

Sextans A (also known as UGCA 205 or DDO 075) From Ott et al. (2012; their Fig. 13), from the VLA ANGST. The high 
resolution (~ 6") H I image shows an elongated (PAhi ~ 135°) disk-like structure, with H I clumps and gaps. The H I morphology 
is approximately perpendicular to the optical axis (PA opt « 45°). The H I is enhanced on either side of the optical disk, roughly 
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positionally coincident with the brightest star formation knots, with opposite extensions to the North and South. The velocity 
gradient direction (PAvg ~ 150°) is slightly misaligned with the H I morphology, but roughly perpendicular to the optical disk. 

UGC 6456 (also known as VII Zw 403) From Thuan, Hibbard & Fevrier (2004; their Fig. 17, 18 and 20), and Felli, Verheijan 
& Fraternali (2014a, b; their Fig. 1, C.12, and Fig. 1, 2), from the VFA. The intermediate resolution (20" x 20", 21.2" x 14.5", 
17.2" x 11" and 15" x 15") H I maps show double-peaked disk emission (PA H i ~ 75°), extended to the Northwest, South and 
West. The H I emission does not align with the main optical emission (PA opt * 115°), but does align with the faint extended 
emission (PA opt ~ 75°). The main H I peak coincides, not with the center of the optical disk, but with the brightest star forming 
regions to the South in the optical image. The velocity field is disturbed, and can be modeled by a combination of rotation and 
radial motions. The kinematical axis (PAy n ~ 90°) is vertical, with the rotation occurring along the H I major axis. 

SBS 1129+576 From Ekta, Pustilnik & Chengalur (2003; their Fig. 1, 3 and 6), with the GMRT. The low resolution (42" x 
40") H I map shows the XMP target, along with its companion galaxy, SBS 1129+577, D ~ 27 kpc to the North, and the dwarf 
galaxy SDSS J113227.68+572142 to the Southeast. The image signals the interaction of the SBS 1129+567/577 pair, by way of 
a H I bridge connecting both galaxies, and a common H I envelope. There is also a small tail of emission to the Northeast and 
Southwest of the XMP. At this resolution, the position angle of the H I emission associated with the XMP target is PA * 75° , 
which is roughly coincident with the direction of the velocity gradient (PAvg ~ 70°). The velocity field is fairly ordered (solid- 
body rotation) in the central region and also in the bridge, consistent with a prograde merger. The high resolution (6" x 14") map 
shows a Northwest-Southeast disk (PAhi ~ 70°), in good agreement with the optical axis (PA opt » 75°). The disk displays small 
extensions to the South, East and West, similar to what is seen at lower resolution. The position of the H I peak encompasses 
both the optical disk center and the bright knot to the Southeast. At high resolution (8" x 7"), the H I peak is resolved into two 
high density regions, each positionally coincident with the optical disk center and the bright star formation knot. Similarly to the 
low resolution image, extended emission to the South, East and West is detected. The somewhat curved tips of the H I disk may 
suggest a warp. The H I position angle on small-scales (PAhi ~ 70°) is roughly in agreement with that on intermediate-large 
scales, with the direction of the velocity gradient and with the optical disk. 

UGCA 292 (also known as NGC 128) From Young et al. (2003; their Fig. 2), with the VFA, and Ott et al. (2012; their 
Fig. 22), from the VFA ANGST. Intermediate resolution (17.7" x 17.4" and 14.2" x 13.9") H I maps show triple-peaked, 
’’boomerang-shaped” H I emission along an almost horizontal axis (PAhi ~ 175°). The internal H I inverted ”V-shape” structure 
mimics what is observed in the optical (PA op , ~ 165°). Indeed, the H I peaks are positionally coincident with the ends and center 
of the optical disk, where bright star formation knots can be found. The vertical gradient direction (PAvg ~ 145°) is misaligned 
by about 20° from the overall H I and optical emission, but is roughly aligned with the Southwest ’’boomerang” optical/H I arm. 
The velocity field displays an asymmetry in the periphery. The high resolution (~ 6") H I spatial distribution and velocity field is 
a scaled-down version of the lower resolution observations. 

GR 8 (also known as DDO 155, UGC 8091 or VIII Zw 222) From Young et al. (2003; their Fig. 3) and Begum & Chengalur 
(2003; their Fig. 2, 3 and 4), with the VFA, Begum et al. (2006; their Fig. 9, 12 and 13,), with the GMRT, and Ott et al. (2012; 
their Fig. 23), from the VFA ANGST. At low and intermediate resolution (30" x 30" and 25" x 25"), the overall H I emission 
is triple-peaked, enclosed in a ’’cross-shaped” halo (PAhi ~ 165°), where the central triple-peak emission is ”V-shaped”. The 
kinematical axis fPA km ~ 80°) is misaligned with both the H I and the optical disk emission (PA opt « 130°). The velocity field 
can be modeled by a combination of circular plus radial motions of the gas, showing distortions in the outskirts. Intermediate 
resolution (18.6" x 18.1" and 14.8" x 14.6") and high resolution (~ 6" and 4" x 3") H I images mimic the triple-peaked ”V- 
shaped” structure (PAhi ~ 165°) seen at lower resolution. These peaks are slightly misaligned relative to the position of the three 
brighter star-forming regions in the optical disk. 

SBS 1415+437 From Felli, Verheijan & Fraternali (2014a, b; their Fig. C.18, and Fig. 1, 2), with the VFA. The intermediate 
resolution (20" x 20") H I map shows a single-peaked Northeast-Southwest (PAhi ~ 115°) disk, extended to the North, and with 
smaller extensions to the South, East and West. The H I peak is positionally coincident, not with the center of the optical disk, 
but with the Southern bright star formation knot that heads the cometary structure, while the cometary tail is mirrored in the 
Northeastern extension. High resolution GO" x 10") H 1 observations show a similar morphology to the intermediate resolution 
image, but where the extended South, East and West emission is now resolved into a web of small tails and filaments. The H I 
emission at intermediate-high resolution scales is roughly aligned with the optical disk (PA opt « 120°) and the kinematical axis 
(PAki n ~ 115°). The velocity field is ordered (rotation) only in the center, with the outskirts displaying lopsidedness. 

Sag Dig From Begum et al. (2006; their Fig. 9), with the GMRT, and Becarri et al. (2014; their Fig. 12 and 13), with the VFA. 
The low resolution (~ 67") H 1 map shows double-peaked quasi-spherical ring-like emission, with an offcenter H I depression 
(AD < 1'). The brightest H I peak roughly coincides with the bright star formation knots to the East, seen on the optical image 
(PA opt ~ 170°). At intermediate resolution (30" x 30"), the ring-like structure becomes more apparent, with multiple H I clumps 
and a large offcenter H I hole. The velocity field displays no sign of ordered rotation. Becarri et al (2014) hypothesized that this 
source could be the result of a merger with another dwarf. 

J2104-0035 From Ekta, Chengalur & Pustilnik (2008; their Fig. 5 and 7), with the GMRT. The low resolution (29" x 25") 
H I map shows a single-peaked inverted ’’triangular” morphology, with several small extensions to the North and a tail to the 
Southeast. The H I peak is not coincident with the bright star formation knot at the Northern edge of the cometary structure, but 
instead is coincident with the Southern edge of the optical disk. The overall H I morphology (PAhi ~ 90°) is slightly misaligned 
from the optical structure (PA opt « 110°), but aligned with the velocity gradient (PAvg ~ 90°). The velocity field displays ordered 
(rotation) in the central part, but a North-South asymmetry is present in the outskirts. The intermediate and high resolution (17" x 
11") map reveals ’’Y-shaped” emission, with extensions to the South, North and Northwest. The overall H I morphology is better 
aligned with the optical structure. In the highest resolution (7" x 5") image, the ”Y-shape” is more pronounced, cradling the 
bright star formation knot in the North. There is extended emission to the South and a clump of emission to the Northwest, 
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mimicking the intermediate resolution image. 



